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Abstract 

We consider two-dimensional system of particles localized on randomly distributed sites of squared lattice 
with anisotropic transfer matrix elements between localized sites. By summing of "diffusion ladder" and 
"cooperon ladder" type vertices we calculated the conductivity for various sites and particles densities. The 
model is relevant to the pjoblem of strong nonmagnetic impurities in superconductors with d x i_ y 2 symmetry 



of the order parameter 

PACS numbers: 74.20-z, 75.10.Hk. 

1 Introduction 

The following tight-binding Hamiltonian is considered: 

H = E *fa " ri)V + (ri)^ + (r,-)p(ri)p(r,-), (1) 

where ip + (ri), ip + (rj) are creating and annihilation operators, p(ri) is the filling number, equal to 1 at the 
localized sites, and to otherwise. The transfer matrix element has a cross-shape configuration 

t(r) = (4,o + Sy, )f(r), (2) 

with 



f(r) = J (^) 7 exp(- 



and a is a lattice constant. We examine the case of a random distribution of impurities on a sites of two- 
dimensional squared lattice. An impurity potential generates a localized state with strongly anisotropic wave 
function. As a result the conductivity is carried out due to hoppings of particles between local states on the 
same vertical or horizontal lines. A similar picture can be realized in 2D g?^^ 2 -wave superconductors, where 
local bound quasiparticle states may arise in the presence of unitary impurities [Q. 

The plan of this article is as follows. In Chapter 2 we consider the case of low impurity density. In Chapter 
3 we calculate the conductivity in the case of high impurity density. In the Conclusion we discuss our results. 

2 Low density 

We consider now the limit of low impurities concentration (c -C 1). In the case of an external electromagnetic 
field we should substitute in ([!]) 



t{r l - Tj) -> t{Yi - rj)exp(ie J A(r,t)dr). 



The electric current is defined as usually from the Hamiltonian (|l|) by varying over a gauge invariant vector- 
potential A 

j a (t) = -ze^(r, - T^ctijV+WVjWiHPj exp(ieA(<)(r 2 - r^)). (3) 
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Where tij — t(r^ — rj). Since we will calculate j{ui) we consider that the potential A depends on time t only. 
Using the equation for the Green function G(t\, t\, r 2 ) 

4-G(t u ri,t 2 ,r 2 ) = -i< T^iil*+(2) > -tf(n - r^fa - h) (4) 
oti 0*1 

we obtain after Fourier transformation in the linear over A approximation 

j a (u) = ?-'52Uj{ri-T j ) a {Ti-T j )i3Afi(w) /'^G(n,r i ,r i )e <na - 

C * — ' / Z7T 

**J*fc.i( r i - r i)a( r fc -r;)^^M y ^G(cj + f7,r J ,r fc )G(fi,r ; ,r i )e jnQ . (5) 

i,j,k,l 

The summation in eqn.(|^) is taken over impurity sites, and a — > +0. In order to evaluate ([5]) in the lowest 
order with respect to the concentration we examine the case of two random arranged sites. The Green function 
is found easy 

G(w, Ti,Ti) = 



(u; + /i) 2 -i 2 / 

G(u,Ti,Tj) = 7— f 1 ' 2 , (6) 



where is the chemical potential. Substituting (^) to (|^) we obtain in the case when sites are arranged on the 
same horizontal chain 

e 2 J 1 
JsM = ^ ( Xl ~ a; 2) 2 ^i,2 ^ + .^ 2 _ 4t2 A a (o>) = Q{u)A x (uj). (7) 

The similar equation can be derived for the case of nonzero temperatures. The result differs only by Fermi 
filling factor. After averaging over impurities sites we obtain for the conductivity a(cu) = iQ(u>)/w 



<j{ui) — —^-c L J x t(x)[np(uj — n) — np{—uj — fj,)]S(u> — 2t(x))dx 
ire 2 2 T XQt(x Q )[np(uj — fx) — np(—u> — /i)] 



't'(zo)| V > 



where np is the Fermi distribution function, 2t{xo) = to. 
Substituting t[x) from (||) we get 

u) exp(-KXo) 
2J = 



(9) 



•'-0 

Tre 2 2 xl[n F {uj - /i) - n F (-o; - /z)j 
8 7 + kxo 

In the limit of low frequency we have following asymptotic behaviors 

1. K = 0, x = (2J/w) 1 /^: 

uj->T,a{u) ocw~ 3/7 , (11) 
ai«T,(TH occj- 3 / 7+1 , (12) 

2. kxo ~ log(2J/w): 

2 J 

w>l>(w) oclog 2 — , (13) 
2 J 

w <C T,a(u) oc wlog 2 — . (14) 

Cl 1 



2 



3 High density. 
3.1 Green function. 

In the case of high density of impurities we will assume that distribution function of impurities may be approx- 
imate by Gausses distribution with dispersion g. That is 



pfo) = C + 5p(Yi) 

< 8p(Ti)5p(Tj) > p = g 2 5ij. 



(15) 



We will assume that concentration c < 1. 

The one-particle Green function for the arbitrary impurity distribution is defined in terms of functional 
integral as usually 



where 



G(E,r,r , ) = (E-t ij p(r i )p(r j ))- 1 r ,= 

J D^D^(r)^(r')exp(iS) 
1 f D$DiPexp(iS) ' 

S = So + Si, 
»S = t5>(iW(r), 

r 

iSi = -i'Y] ^(ri)i(ri - r 2 )p(ri)p(r 2 )^(r 2 ). 



(16) 

(17) 
(18) 
(19) 



Introducing an additional integration over new fields \ > X m order to eliminate the second order terms pp, we 
get 



e lSl = y^x^Xexpj^^x^i)^ 1 ^! -r 2 ) X (r 2 ) 
+c^(x(r)V(r)+^(r) X (r)) 

r 

+ ]T <5p(r)( X (r)V(r) + ^(r) X (r))}/Z 

r 

Z= DxD X exp{-i^2 x^i^Hri -r 2 )x(r 2 )} 

J r ir2 
k 

e(k) = ^t(r)e^ = -Jln[( K 2 + 4 sin 2 ( K 2 + 4 sin 2 (^))] 

The Green function in terms of new two component field ip 



where 



(20) 

(21) 
(22) 
(23) 



^i0) 

^i( r ) 



X(r) 
^(r) 



reads 



(24) 



G(ri, r 2 ) = -i < <p(ri) ® £(r 2 ) > 



(25) 
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where ip(r) = (<£>i(r), y 2 (r)) and angle brackets are denned as 

/ D$Dil)DxDx{- ■ .)e iS °ff 



<...>- 



with 



/ D^Dil)DxDxe iS eft 
iE c + 8p(r) 



c+ 8p(r) -it 1 (r 1 - r 2 ) 
The equation for the Green function after averaging over impurities in the Born approximation reads 

G(k) = G°(k) + G°(k)E(k)G(k) 

with the bare Green function 

= [ \ :f- 1(k) 

and the self-energy S obtained by summing of diagrams without intersections 

dki 



S(k) = g 2 a 2 



The solution of Eqns. (E8h, (BOh is 



E(k) = 



Q -iP 



a x Gl}n)u x . 



-iP R 



G(k) = 



1 



1 + e(k)P)(P -Q)-(c+ P) 2 e(k) 



l + e(k)P -i(c + P)e(k) 
-<(c + P)e(k) -e(k)(P-Q) 



where 



p 
p 



2 2 

9 a 



dk 



(P - Q)e(k) 



2 2 

"5 a 



2 2 

3 a 



(2tt) 2 (1 + e(k)P)(P -Q)-( c + P) 2 £ (k) 
dk 1 - e(k)R 

(2tt) 2 (1 + e(k)P)(P - Q) - (c + P) 2 e(k) 

dk (e + P)e(k) 

(2^) 2 (1 + e(k)P)(P - Q) - (c+ P) 2 e(k) - 



In the limit of low dispersion g 2 <C 1 we obtain 



Q 



p 





2 


R,A _ 


9 

tr 7 
T ^2 


R,A 


±41 

P 2 


= 2?rg 2 a 2 


P 2 (E\ 



(26) 
(27) 

(28) 
(29) 

(30) 

(31) 
(32) 



Where ^o( £ ) is density of states of the pure model (c=l, g = 0): 

r jit 

V{£) =] (2^ <S(£_£(k)) - 
Taking into account that Q <C P, PP C 1, P <C c, we find for the Green function in the limit j 2 <1 



G R ' A (k) 



1 — ice(k) 
P-c 2 e(k)±i 7 /2 V -tce(k) -e(k)P 



(33) 



(34) 



(35) 



(36) 



4 



3.2 Drude formula 

The conductivity in our case is defined as in Ch.I in terms of four-particle correlation function 

? 2 f dki dk 2 



a E(u) = |^ J ,^f 2 v a (ki)v ot (k 2 )KEuO<-i, k 2 ; k 2 , ki) 



where E is taken at the Fermi level, 



ifE W (ki,k 2 ;k 2 ,ki) = 

1 £ e ^(x-y) e ,k 2 (z^t) (pxPyPzptG ^ (y5Z ^ + w/2)G A (tjX ^_ Lj/2))p 



x,y,z,t 

and 



Substituting the solution (36) to (mm we find in the lowest approximation 



p 2 r jk 2 6 / T 

M^) = ^ y ^^(kJGft (k, £ + W /2)G&(k, £ - W /2) = — 

where 



c 6 fjy A(E)_ 
B(E) 



A(E) 



[ %^(k), S(E) = / 



v(k) — ^Jv^Ck) + u^(k),and dl^ - element of the length of the Fermi surface. 
The conductivity can be expressed in terms of a particle density defined as 



/7r / a r7t I a 

dk x / dk y 9(E/c 2 - e(k)). 
-7r/a J— n/a 



dk x / 

T/a J—n/a 

We obtain in limits of low and high densities following results: 



7 n , for n <. 1, 



2 c 6 

P~327rlog 2 (2)' 

^4^i5Fw( 1 " n o)' for (l-«o)«l. 



<J 



the asymptotic behavior in intermediate region < no < 1 is 

P 2 6 1 



9 2 (l-no^log^) 
with maximum value 



maa; K 8/3 fog(K) 

reached at 1 — no ~ k 4 / 3 . 

3.3 The absence of weak localization. 

To go beyond the quasiclassical approximation we include contributions to the conductivity from 
ladder" and " crossed-ladder" or " cooperon" vertices || . The addition term is 

5<j e (cj) = ^-J ( 2^ 2 |^ «a(k 1 )^(k 2 )Gf g (k 1 )Gf l (k 2 )Gf fe (k 2 )G^(k 1 )A-ac; M (ki,k 2 ) 
where for "diffusion" - vertex contribution we get the equation (see Fig. la): 

fr(-°) — ^2x x ,2 f d\ _ R /,\t^(D) n A ft\„x 
n ac;bd ~ 9 a ac a bd+ 9 J x 2 1 1 Cl c l c ; bb l 61 dl ^ ' dl d ' 



The solution of this equation does not depend on ki and k2 . Therefore the contribution of a " diffusion" - vertex 
to the conductivity is equal to zero because 



a) 



/ (^2^( k l) G fa( k lX( k l)=0- 









b) 











(48) 




Fig. 1. Diffusion a) and crossed-ladder b) vertices . 

Now we consider the "cooperon" -vertex contribution. The vertex K' D ' obeys the equation (see Fig. fb ) 

K { £,M) = 9 2 <?> X bd + 9 2 J J^<e°lA (l)Gf bl (q - I)<2w(q), ^ 
where q = ki + k2 . 

If the solution of this equation has some singularity behavior as function of q (for example " diffusion pole" ) we 
may hope to obtain the contribution to conductivity. Now we prove that "cooperon" -vertex has not singularity. 
We search the solution of in the form 



js\ c ) _ n 2 \ " n v a 

^■ac-bd — 9 2-^ a,c a bd A ^i>i 



(50) 



where \x^v = 0,x,y,z, and A satisfies the equation 



a(3 



with 



Aa p = 9_ [ J}_ 

» v 4 J (27T) 2 



[Sp((j x G R ^a a )\ [Sp{a x G R a v a )] 



The matrix A can be rewritten as 
with 

and 



A af) = -A a A p 



i E/c 
E/c -iE 2 /c 2 



(51) 

(52) 

(53) 
(54) 

(55) 



G 



So the components of matrix A are 



A = 



/ 1 isinri(0) cosh(6>) 

i sinh(6>) 1 

cosh((9) 1 

y icosh(6>) sinh(6») 





-i cosh(6') 
- sinh(0) 
1 



where 



c/E = e e 

The solution of (p2) is found in terms of matrices U, B: 



where 



1 



Bui, 



The eigenvalues z v are following 



1 + z p z v 
det(A - zl) = z 2 (z - 2) 2 = 



^o = 

zi = 

z 2 = 2 

z 3 = 2 



(56) 

(57) 

(58) 
(59) 
(60) 
(61) 



(62) 



As a result we have that for q = the function B and consequently K has no singularity. Therefore we do not 
obtain the phenomenon of "weak localization" in our 2D model due to specific type of a disorder. 



4 Conclusions 

We have investigated the two-dimensional model with a new type of disorder due to a random distribution of 
local states with strongly anisotropic overlaps of wave functions. The conductivity of this system was calculated 
in limits of low and high densities of local states. We have shown that considered type of disorder does not lead 
to the weak localization phenomenon, as usually in two-dimensional case [0, due to the absence of logarithmic 
divergence from the integration over the diffusion pole. 

5 Acknowledgments 

S.M. thanks A. Balatsky for a discussion. This work was supported by the Russia Fund of Fundamental Research 
under grant No 960217791. 

References 

[1] L. P. Gor'kov, A. I. Larkin, D. E. Khmelnitskii, Pis'ma Zh.Eks.Teor.Fiz. 30, 248 (1979). 
[2] P. A. Lee, T. V. Ramakrishnan, Rev.Mod. Phys. 57 , 287 (1985). 

[3] E. Abrahams, P.W.Anderson, and T.V. Ramakrishnan, Philos. Mag. B 42 ,827 (1980). 

[4] A. V. Balatsky, M. T. Salkola, Phys.Rev.Lett. 76, 2386 (1996). 

[5] 



7 



